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Abstract 

We consider solutions where the surface tension of the solvent 7 S is 
smaller than the surface tension of the polymer y . In an evaporating 
film, a plume of solvent rich fluid, then induces a local depression in 
surface tension, and the surface forces tend to strengthen the plume. We 
give an estimate (at the level of scaling laws) for the minimum thickness 
e* required to obtain this instability. We predict that a) e* is a decreasing 
function of the solvent vapor pressure p v (e* ~ Pv 1//2 ) b) e* should be 
very small (< 1 micron) provided that the initial solution is rather dilute, 
c) the overall evaporation time for the film should be much longer than the 
growth time of the instability. The instability should lead to distortions 
of the free surface and may be optically observable. It should dominate 
over the classical Benard-Marangoni instability induced by cooling. 

PACS numbers: 68.60Bs; 68.45Da; 6810 Jy. 

Shortened version of the title: 

Film evaporation 

1 Introduction 

Spin casting of polymer films from solutions is an important practical process. 
The intrinsic state of these films after casting raises a number of questions - 
especially if they are thin In most cases, we are dealing with materials such 
as polystyrene, or polymethyl metacrylate, which are glassy in the final state: 
these glassy features complicate the process enormously. 

In the present pages, we consider a simpler case: a) the dry polymer is 
assumed to remain fluid at room temperature (for instance, it could be a silicone 
oil) b) the polymer is assumed not to adsorb on the free surface of the solutions 
(this corresponds to 7 p >7 s ). We shall see that certain convective instabilities 
should occur during evaporation: one driven by concentration effects, and one 
driven by thermal effects (Benard Marangoni) . We find that the concentration 
effects should usually be more important. 

In section 2, we give a simple description of the classical evaporation process 
in a mixed film, associated with a diffusion flow of the solvent. Many more 
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precise discussions of this process exist in the literature ||, but the present 
crude picture is enough to estimate the concentration drop ip d — ip u between the 
bottom plate (ip d ) and the upper free surface (ip u ). This drop is the essential 
control parameter for the onset of a convective instability (section 3). 

This instability requires low viscosities: indeed the films are initially dilute 
and of low viscosity. We discuss the dilute limit in section 4. 

Ultimately, we discuss another possible source of instability related to ther- 
mal effects: evaporation cools the outer surface. This leads to a classical Benard 
Marangoni instability j|. We show in the appendix that, (when j p > 7 S ), con- 
centration effects should dominate over thermal effects. 

All our discussion is restricted to the level of scaling laws. This is not unrea- 
sonable: even without instabilities, the concentration profiles in the evaporation 
problem are quite complex, and depend on many details: an exact hydrodynamic 
calculation of the thresholds in this situation would be purely numerical, and 
not very informative. 



2 Steady evaporation profiles 

The volume fraction of solvent ip{ z ) decreases when we move from the bottom 
plate (z = 0, ip = tpd) to the free surface (z = e(t), ip = ip u ). At the free surface, 
we assume instant equilibrium: the partial gas pressure of the solvent is equal to 
the equilibrium gas pressure p e (tp u ), corresponding to ip u - I n practice, we shall 
be concerned mainly with dilute solutions and replace (very roughly) p e (ip u ) by 
the vapor pressure of pure solvent p v (T). Above the free surface, the solvent 
molecules diffuse in air, with a certain diffusion coefficient D a . In practice, the 
air transport involves both diffusion and convection, dependent on various noise 
sources in the experimental room. As is often done, we assume that all this 
can be described as diffusion through a boundary layer of fixed thickness £ (~ 
1mm). The number density of solvent molecules in the gas 7g(z) thus drops 
linearly from the value associated to p v : 



Jc(z = e) = ^ (1) 
to at the upper end (z = e + £). 

The diffusion current in the film (number of molecules per unit area and per 
second) is: 



W = D coop ^ a 3 ^ (2) 
are 

where D coop is the cooperative diffusion coefficient of the solution ||, and 
we ignore the tp dependence of D coop . The same flux is found just above the 
free surface: 



w = Da la(e)-Ja(e + t) =Da Ja^ (3) 

Our scaling estimate for the diffusion constant of solvent in air is D a i r ~vt^A, 
where v t h is a thermal velocity for a solvent molecular, and A is a mean free 
path, inversely proportional to the air pressure p a . This gives ultimately: 
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D air £* -f — 4 



where a is the size of a solvent molecule. 
Combining eqs(|l[]|), we arrive at: 



^d-^u~^ (5) 

^ coop ya 

Typically, taking e = £ = 1mm, av t h = 10~ 3 cm 2 /sec, D coop = 

10~ 6 cm 2 / sec, and p v = 0.1p a : we obtain ^ d — -0 U ~ 0.1. Eq. (||) is important, 
because it defines the driving force for the convective instability to be discussed 
in the next section. 



3 Concentration 

The principle is shown on fig.|]. During the unperturbed evaporation process, 
we saw that there is a significant difference between the solvent fractions ip d 
and ip u (down and up) at the supporting surface and at the free surface. We 
now assume that a roll instability is superposed on the diffusion flux. We can 
visualise the system as a set of rolls of size e, or as a set of plumes carrying extra 
solvent from down to up. The result is a slight difference in solvent concentration 
between point A (4>a) an d point B (4>b)- We estimate this difference by a balance 
between convection and diffusion: 



D coo / A ^ = Vfa - (6) 

e 

where V is the convective velocity. The difference ip A — ip B induces a gradient 
of surface tension: 



V7= = Wa-tPb) ( 7 ) 

e e 

where r y' s =d'y / dip is assumed negative: the polymer has a higher surface 
tension than the solvent. 

Finally, we estimate V by a balance between surface forces and viscous 
stresses: 



V 7 ~ n — (8) 

e 

where rj is the solution viscosity. 

Combining eqs 0, |^) one arrives at a minimal thickness: 



* ^ yDcoop , q s 



3 



Figure 1: 



For rough estimates, we can write in the semi dilute regime: 

kT 

Dcoop = —z (10) 

where kT is the thermal energy, n s the solvent viscosity, and £ the mesh size 
in the solution || . 
The result is: 



kT ri 

(11) 



\l's\^Vs(^d-^u) 

We can now combine eq. |ll| and eq. |^, arriving at: 



(e*) 2 ~ J 1 £n Dcoop Pa 
\l'\S,Vs av th Pv 



4 Discussion 



a) Role of the solvent vapor pressure p v : we see on eq. (12) that high 
values of p v lead to low values of e*: the instability (which occurs for e > e*) is 
favored. 

b) Role of the solution viscosity 77: if we were dealing with concentrated 
polymer hims, r]/rj s would be very high, and e* would be prohibitively large. For 



mm, 



instance, if we choose |7'| = 10mJ/m 2 , £ = lOnm, i]/r] s = 10 6 , 1=1 
and D coop /vthd — 10~ 3 , we arrive at e* ~ 0.2mm. However, we must not forget 
that our films always start from a rather dilute state, with r]/r) s ~ 1. The 
instability will grow during the early stages of evaporation. If we now switch 
to rj/rjs ~ 1, we are led to e* = 200 nanometers. Thus (even if there are large 
numerical prefactors in eq. (|l2])) we expect to find the instability in films around 
1 micron. 
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c) To substantiate this prediction, we must show that the evaporation 
time T ev is longer than the growth time (r g ) of the instability. We estimate r cv 
as follows. 

The rate of thinning is given by the flux W of eq. (0) : 



I = - 3w < 13 > 



Using also eq. (g), this gives an evaporation rate: 
1 1 de D, 



■'coop 



T ev e dt e 
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(14) 



Lest us now consider the growth time of the instability t g . For simplicity, we 
shall focus our attention here on rather thin films, where inertial effects can be 
neglected. We can then rewrite a time dependent equation for the modulation 
ip m = ip A — ip B in the form: 

^ = - Dcoopq ^ m + q V(^ d -i; u ) (15) 

where q is the horizontal wave vector of the modulation (q ~ e _1 ). We also 
have, as in eq. (@): 



m V = q\i\^ m (16) 



At threshold, the two terms in eq. ( |15| ) balance each other. We now consider 
situations definitively above threshold, where the convection term in (15) is 
dominant. Using eq. [l6|, this leads to: 

l IV I 

-e^m^-VO (17) 

Tg -qe 



Comparing eqs ( p7| ) and (14) we find: 

T 9 ?]Dcoop 



(18) 



Taking a dilute system, with j^'l = IQmJ/m 2 , e = lmicron, £ — 10nm, 
we find T ev /r g ~ 2.10 4 . Thus, the instability should grow easily before complete 
evaporation -during a time interval r g << r ev , where the viscosity is low and 
e* is small. 

Another estimate for r g , can be written down, when inertia is dominant (and 



eq. 16 is modified). But the conclusion remains the same for most practical film 
thicknesses. 

d) There are a number of details which are not included in our dis- 
cussion: for instance, the dependence of the surface tension ^/{tp) on solvent 
concentration in the dilute limit (ip — > 1) has not been considered. In fact, we 
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know that when the polymer surface tension is higher than the solvent surface 
tension, there is a depletion layer of thickness £ near the free surface, and the 
surface tension has the following scaling form || : 

7W - 7 (i)-^-§(i-</0 3/2 (19) 

4 a 

Thus, the parameter |7'| of eq. (Q) is really ip dependent: 

| 7 '|~^(WJ 1/2 (20) 

But, for reasonable values of ip u (e.g. ip u = 0.9) this brings in only minor 
corrections on e* -proportional to {l-ip u ) 1 l i . 



5 Conclusions 



D t — = V(T d -T u ) 



(Al) 

The difference T d - T u is 
The vertical solvent flux 



All our discussion ignored numerical coefficients, and we know how they can be 
important (e.g. for the onset of thermal convection in the Rayleigh problem). 
But some conclusions do emerge: 

1) The convective instability due to concentration gradients in the film should 
show up for non glassy polymer films (of thickness e > e*). 

2) ft takes place early: when the film is still dilute. 

3) It will induce a certain surface roughness. To compute the amplitude of 
the surface undulations resulting from the rolls is a delicate task and will not 
be attempted here. 

Acknowledgments: we benefited from discussions with G. Reiter, S. Kumar 
and F. Brochard-Wyart. 

6 Appendix 

We now ignore all concentrations effects, and focus on fig. |[ Here, we have 
a thermal plume, terminating at point A, and we estimate the temperature 
difference by a balance of heat fluxes: 



e 

D t is the thermal diffusion coefficient of the liquid, 
due to evaporation and cooling at the upper surface, 
towards the air is: 



J s — D COO p (A2) 
e 

This corresponds to an outward energy flux: 
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Figure 2: 



Je = U v J s (A3) 

where U v is the heat of evaporation (per unit volume of liquid) . 
Balancing Je against a vertical heat flux (from the support to the outer 
surface), we obtain: 

C p D t ^—^ = J E (A4) 

e 

Here, C p is the specific heat per unit volume (assumed constant), and D t a 
thermal diffusion coefficient in the liquid. 

Finally, we have a stress balance similar to eq. |8[ giving: 

r]V £* -~f' T (T A - T B ) (A5) 

where: 

It = % < (A6) 
Combining eqs (A2-A5), we arrive at a threshold thickness: 



D 2 t C p 1 



coop 

u 

We can now compare the two processes, however: 



(A7) 



e 



is c p f D 



e* i T U v \D. 



■'coop 



(A8) 



The factors 7gC p /7^J7 v combine to give a constant of order unity, and 
the major features are related to the diffusion coefficients. Typical values are 
D COO p = fCT 6 cm 2 / sec, D t = fCT 3 cm 2 /sec, and we thus expect e** >> e*: for 
films with 7 p > 7 S , thermal processes should be dominated by concentration 
processes. 
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